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ABSTRACT 

The  consequences  of  assuming  that  density  is  conserved  in  the  problem  of  determining 
absolute  velocities  are  investigated.  Two  questions  are  considered:  (i)  the  constraints  that 
the  density  must  satisfy  to  be  compatible  with  the  assumed  geostrophic  and  hydrostatic 
dynamics  and  (ii)  whether  and  to  what  extent  the  indeterminacy  in  this  dynamics  is 
removed  by  this  additional  assumption. 


DTIC 

ELECTEI 


NOV  2  3  1987 


DISTRIBimON  STATEMENT  A 


Approved  fc>r  priJjlic  reloraw: 

IX.itnb ‘io  -.  '  i !• 


1  1  6 


1.  Introduction 


'With  the  dawn  of  ocean  forecasting  approaching,  the  need  to  know  velocity  fields  for  nse 
as  initial  conditions  in  a  prediction  scheme  is  keenly  felt.  Bat  direct  velocity  measurements 
are  still  difhcalt  to  make,  and  so  we  are  forced  to  rely  on  the  classical  procedure  which 
consists  of  inferring  the  velocity  fields  from  temperature  and  salinity  measurements.  In 
this  procedure,  we-arV  given  the  density  field  .pf»r^,^)^'amHN^mast  find  velocity  fields 
v(x,y,z),v(x,y,2),u;(z,y,z)  and  a  pressure  field  p(x,y,  z)  such  that 


fv  =  PoVi 

fu  =  -Pc^Pj. 

0  =  -p"^p,  -  p 
li,  +  Vj  +  u),  =  0 


(11°) 

(1.16) 

(1.1c) 

(lid) 


The  nutations  in  the  above  equations  are  standard:  x,y,z  are  cartesian  coordinates;  /  is 
the  Coriolis  parameter  and  g  is  the  gravitational  acceleration:  po  is  an  average  value  of  the 
density. 


As  is  very  well  known,  the  integration  of  the  hydrostatic  equation  from  the  ocean’s  bottom 
at  z  =  /i  to  a  depth  z,  viz. 


p=  -9  J  -t-  n(x,  y) 


introduces  an  unknown  function  of  x,y  denoted  here  by  11.  This  ‘barotropic’  component 
of  the  pressure  cannot  be  determined  from  the  dynamics  (l.la,b,c.d). 


The  indeterminacy  in  the  pressure  entails  a  similar  indeterminacy  in  the  velocity  fields. 
Indeed,  if  we  define 

P{x,y,z)  = -g  J  ^p{x,y,z')dz'  (1.3) 


u  =  -(oof)-^{p,-^n,} 
V=  (pofrUPr  +  Ur) 


(1.4a) 

(1.4h) 


The  vertical  velocity,  which  is  obtained  by  integrating  the  continuity  equation  in  the  ver¬ 
tical, 

w  =  w{x,y,-h)-  J  {u^  +  v,)dz'  (1.5) 

abo  contains  the  same  indeterminacy.  The  search  for  the  famous  'level  of  no  motion’  is 
intimately  related  to  the  determination  of  n(x,y). 


Following  in  the  footsteps  of  Worthington  (1976),  Wunsch  (1978),  Fiadeiro  I:  Veronis 

(19821  and  others,  we  want  to  examine  whether  the  addition  of  extra  information  in  the*'^ - - 
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form  of  th«  conservation  of  the  density  field,  vis. 


up,  +  vp,  +  wp,  =  0  (1.6) 

can  determine  n(x,y). 

Mathematically,  it  is  not  clear  whether  Equation  (1.6)  is  compatible  with  the  dynamics 
(1.1).  We  shall  examine  this  question  first  on  the  /—plane  and  then  on  the  plane. 
We  shall  show  that  in  both  cases  the  density  field  must  satisfy  certain  constraints.  If 
these  constraints  are  satisfied,  then  for  the  ^-plane  case,  but  not  for  the  /-plane,  II  can 
be  determined.  Our  analysis  is  carried  out  in  the  simplest  case  of  a  flat  bottom  ocean. 
However,  it  can  be  extended  to  the  case  where  bottom  topography  exists  as  well  as  to 
cases  where  quantities  other  than  density  are  conserved. 

2.  A  mathematical  result 

An  essential  step  in  our  analysis  consists  in  finding  a  z— independent  solution  of  an  equation 
of  the  form  ' 

A(i,y,z)ni  +  5(i,y,z)n,  =  1  (2.1) 

which  holds  in  a  cylindical  basin  D  =  {(*.y,z)  :  (x.y)  £  D,—h  <  z  <  0},  where  D  is  the 
horizontal  cross-section.  If  A  and  B  themselves  are  independent  of  z,  then  the  existence  of 
such  solutions  is  not  too  surprising.  Of  greater  interest  is  the  case  in  which  such  solutions 
exist  eventhough  A  and  B  are  z  dependent,  i.e  when 

Al->-Bl^0  in  D.  (2.2) 

Relabeling  coordinates  if  need  be,  we  interpret  (2.2)  to  mean  that 

B,^0  in  D  (2.2') 

.  With  the  above  conditions  we  have  the  following  result. 

Theorem:  There  exists  a  z-independent  solution  n(x,y)  of  (2.1)  if  and  only  if  for  all 
(x.  y.  z)  (r  D 


(i)  A{x,y,z)  =  koix,y)B{x,y,z)  +  lo{x,y)  (2.3a) 

(u)  /o(x,y):^0  (2.36) 

(iu)  =  0  (2.3c) 

Furthermore,  this  solution  n(x,y)  satisfies 

n.-  =  /„-'(iCf)  (2  4o) 

n,  =  (2.«) 


Proof:  Suppose  (2.3a,b,c)  hold  in  D.  We  must  show  that  a  z-independent  solution  n(z,  y) 
of  (2.1)  exists,  and  that  it  satisfies  (2.4a, b). 

In  view  of  (2.3b),  we  can  write  (2.3a)  thus: 

=  l  (2-5) 

Now,  (2.3c)  implies  the  existence  of  a  function,  say  ll(z,y),  such  that  (2.4a, b)  hold.  Sub¬ 
stituting  /q  '  and  -ho/Q  ^  by  their  expressions  in  terms  of  IT  in  (2.4),  we  conclude  that  this 
function  11  is  a  solution  of  (2.1). 

Conversely,  suppose  that  11  is  a  z-independent  solution  of  (2.1)  and  that  (2.2*)  holds.  Then 
we  can  differentiate  (2.1)  with  respect  to  z  and  write 

=  0.  (2.6) 

W'e  note  at  this  stage  that  A^B  -  B^A  ^  Q  in  D.  Indeed,  if  this  were  the  case,  then 

B.iATlr  +  fillj,)  =  Bf. 

which  is  obtained  by  multiplying  (2.1)  by  B^,  could  be  written  as 

B{AtBg  •+•  .Bjlli,)  =  Bf. 

But  as  a  result  of  (2.6),  we  reach  a  contradiction.  Therefore 

A^B-B,A^0  (2.7) 


and  as  a  result,  we  can  solve  (2.1)  and  (2.6)  for  I!,  and  11,  to  find 


n,= 


AB,  -  BA, 

A, 


AB,  -  BA, 


(2.8a) 


(2.8b) 


This  step  is  not  nn^^e  that  taken  by  Stommel  &  Schott  (1977)  and  Needier  (1985)  in 
their  study  of  this  problem. 

W'e  must  next  insure  that  these  expressions  for  11  are  indeed  z-  independent.  By  considering 
their  ratio,  we  immediatedly  see  that 


A,  -  ko(x,y)B, 


and  after  integration  we  arrive  at  the  condition  (2.3a)  between  A  and  B.  Substituting  this 
condition  in  (2.8a.b),  we  arrive  at  (2.4a.b).  Condition  (2.3c)  arises  from  forming  Ilxy  in 


two  dilTerent  ways.  Finally,  the  condition  (2.7)  translates  into  (2.3b).  We  can  easily  verify 
that  farther  z  differentiations  do  not  introduce  other  conditions.  This  completes  the  proof. 


In  closing  this  section,  we  give  an  example  of  a  partial  differential  equation  of  the  form 
(2.1)  in  which  the  coefficients  A  and  B  satisfy  all  the  conditions  of  the  theorem.  The 
example  is: 

fS(*,v,2)  +  -777^^ - rin*  B(i,v,z)n-  =  1 

j>{x-yr 

where  B  is  an  analytic  function  of  2.  In  this  case, 

^0  =  1 

/'(®  -  y) 

and 

n  =  /(*  -  y) 

3.  The  /-plane 

Wc  consider  in  this  section  the  case  of  the  /-plane,  simply  as  a  foil  for  the  /5 -  plane  which 
is  discussed  next.  In  this  case,  the  flow  is  horisontally  non-divergent.  This,  together  with 
the  fact  that  the  bottom  is  flat,  implies  that 


ID  =  0 

Thus,  the  added  conservation  equation  reads 

up*  -I-  vpy  =  0 


(3.1) 

(3.2) 


Making  use  of  the  hydrostatic  and  geostrophic  balances,  this  equation  can  be  written  as 
follows; 

~P»P**  'h  P*Pjr*  —  0  (3.3) 

and  since  pi  -t-  is  not  identicrdly  sero  in  D  lest  the  ocean  would  be  at  rest,  this  implies 
that 


’■  (s)  = 


or  equivalently 


P*  =  -P,  tanQo(x,y) 


v(x,l/,z)  =  u(x,y,z)Unao(x,y) 


Thus,  at  a  given  latitude  and  longitude,  the  velocity,  which  is  purely  horizontal,  has  the 
same  direction  for  all  depths!  This  result,  in  turn,  places  some  constraint  on  the  data. 
Indeed,  substituting  (3.5)  in  (3.2)  shows  that 


This  constraint  is  typical  of  the  constraints  the  conserved  quantity  must  satisfy  for  the 
data  to  be  compatible  with  the  assumed  dynamics.  What  this  means  in  practice,  is  that 
at  each  station,  2— independent  ao  are  generated  by  passing  the  best  fitting  line  passing 
through  points  with  absiscas  and  p,  (or  whatever  the  corresponding  quantities  are  for 
the  case  at  hand).  Figure  i  shows  an  attempt  to  determine  one  such  angle  ao  from  data 
from  the  REX  (Regional  Energetics  Experiment)  taken  on  August  13,  14  and  15  1985  at 
a  station  in  the  North  Atlantic  located  at  latitude  37^  N  and  longitude  64^  E.  The 
density  was  obtained  from  AXBT  temperatures  measured  every  two  meters  up  to  a  depth 
of  743  meters  merged  with  climatological  salinity  data.  Figures  2(a),  (b)  and  (c)*show  the 
same  data  segregated  by  depth. 

As  we  shall  see,  a  similar  constraint  arises  on  the  /?-plane. Having  satisfied  this  constraint, 
are  we  now  capable/io  determine  11?  The  answer  is  unfortunately  no.  Indeed,  from  the 
definition  (1.3)  of  P.  it  follows  that: 

P,  +  P,  tan  ao  =  0  (3.7) 

This  means  that  the  baroclinic  components  of  the  flow  field  satisfy  the  conservation  equa¬ 
tion  by  themselves.  Therefore,  the  need  for  a  barotropic  correction  is  not  imperative.  From 
f3.4)  and  (3.7)  it  follows  that 

IIx  +  Ily  tan  ao  =  0  (3.8) 

If  we  were  forced  to  proceed  further,  then  we  would  need  to  know  the  distribution  of 
sources  and  sinks  at  the  edge  of  out  basin  in  order  to  solve  this  equation  for  H. 

4.  The  plane 

Because  of  the  /S-effect.  the  horizontal  mass  divergence  b  —3vff(\  and  therefore  the  vertical 
velocity  as  given  by  (1.5)  is 


w{x,y,  z) 


i  t\  1  $ 

v(x,y,z  )az 


(4.1) 


or  in  terms  of  P  and  11 


=  {p/fopo)i  f 


J-h 


+  [z  -t-  ftjll,/ 


We  shall  find  it  convenient  at  times  to  write 

W  =  {(3if^Po)  f  j^dz' 
The  conservation  equation  for  p  now  becomes 


{p,  +  (0\z  +  h]/fo)pt}U,  -  =  -p,P,  +  p,P,  -  {pofo)ptW 


(4.2) 


(4.3) 


(4.4) 


*  We  are  indebted  to  Dr.  Ted  Bennett  for  providing  ns  with  these  figures. 


The  question  once  again  is  whether  the  right  hand  side  of  (4.4)  vanishes,  i.e.  whether  the 
baroclinir  fields  identically  satisfy  the  conservation  equation.  We  recall  that  for  0  =  0  the 
answer  to  this  question  was  always  yes.  If  ertch  side  of  (4.4)  is  separately  tero,  then  the 
density  field  which  makes  up  the  data  is  constrained  to  satisfy  two  distinct  equations.  The 
equation  stemming  from  the  baroclinic  part  is 

-  (po/o)p*W^  =  0  (4-5) 

whereas  that  from  the  barotropic  part  yields 

P*  +  {P»  +  (/5(^  +  ^]//o)Pi}  =  0  (4.6) 

where  ao  is,  as  previously,  solely  a  function  of  x  and  y.  Are  there  isopycnal  surfaces  which 
satisfy  both  of  the  above  equations? 

We  shall  show  that  this  question  must  be  answered  negatively.  Integrating  (4.6)  over  z 
and  using  the  dehnition  (1.3)  of  P.  we  deduce  that 

i-Unao{P, -t  yP)  =  0  (4.7) 

JO  JO 

Next,  we  eliminate  p,  and  Py  from  (4.5)  using  (4.6)  and  (4.7).  The  result  is: 

fs  f  <‘p} = *’}  (<■«) 

If  we  differentiate  this  equation  with  respect  to  z,  we  see  that 

PjQ^Pxt{P  -r  pg  f  dz'}  =  0fo'^p„{-Px  f  dz  ^  f  Ptdz'}  (4.9) 

J-h  j-h  J-h 

The  above  two  equations  imply  that 


PtPxM,  —  pxptt 


i.e. 

P*  =  pxtan7o(a:,!/)  (4.10) 

V»e  have  thus  reduced  the  problem  to  finding  whether  there  exists  a  function  p(x,y.z) 
which  satisfies  the  system  of  equations 


E\\p]  =  tan-voPr  -  Pi  =  0 

E'i\p\  =  Px  +  tauQop,  +  (0\z  +  /il//o)tanQop,  =  0 

At  this  stage  we  appeal  to  the  theory  of  systems  of  first  order  partial  differential  equations 
(see  e.g.  Smirnov  1964,  p.34.).  The  first  step  in  the  theory  consists  in  attempting  to  derive 
further  independent  equations  for  p.  Such  additionrd  equations  are  obtained  by  forming 
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PoisBOD  brackets.  Thus,  denoting  temporarily  z,y,2  by  zi,  Z3,  zj  and  p.,  p,,  pt  by  pi,  pj, 
P3,  the  theory  leads  ns  to  consider 


3 

E,[p]  =  53{(«f:,/a*.)(«£,/ap,)  -  (dE,/dx,)(dE,fdp,)}  =  0 
1-1 


l.e. 


-[(^-ro/^*)  +  (57o/^v)tanao)(cos'ro)*^]Pr  +  [(^ao/5y)tan7o(co»cto)”^]Py 

-t-  [/3(2  +  h)/"^(dao/dy)tan7o(co8ao)~^  —  3f~^  tanaolpi  =  0  (4.11) 

This  equation  is  independent  of  the  other  two.  Having  obtained  three  independent,  homo¬ 
geneous.  first  order  partial  differential  equation  for  p  which  is  a  function  of  three  indepen¬ 
dent  variables,  the  theory  tells  us  that  p  =  constant  is  the  only  possible  solution. 

Thus,  on  the  /i-plane  the  baroclinic  velocity  fields  cannot  satisfy  the  density  conservation. 
Therefore,  we  can  write  (4.4)  as 


where 


and 


-4!U  BUy  =  1 

(4.12) 

,  ^*+p,f^(2  +  h)//o] 

P.P,  -  p,P.  -  P.w  ’ 

u 

(4.13a) 

B  - 

PxPf  -  P,Pt  -  PtW 

^  fH 

(4.136) 

If  we  assume  that  A  and  B  are  functions  of  depth,  then  the  theorem  yields  the  following 
constraint  on  the  data 

A  =  ko{x,y)B  +  lo(x,y) 

and  the  usual  expressions  for  the  barotropic  components,  namely 

n.=/o'(z,y) 

n,  =  -h)(x,y)/lo{x,y) 

On  the  /0-plane,  the  density  conservation  is  stringent  enough  to  determine  the  barotropic 
field  completely. 

5.  Concluding  remarks 

We  have  attempted  to  elucidate  the  question  of  how  a  depth  dependent  conserved  field 
ran  remove  the  depth  independent  indeterminacy  of  the  velocity  fields.  We  have  seen  that 
this  attempt  leads  naturally  to  a  characterization  of  the  class  of  acceptable  density  fields. 
In  turn,  this  suggests  a  very  natural  variational  procedure,  namely  that  the  “distance” 


between  the  data  and  this  set  of  acceptable  densities  be  mimimiied.  For  the  plane, 
there  is  a  nniqne  Telocity  field  associated  with  this  projection.  Mr.  T.  Bennett  of  NORDA 
is  currently  implementing  this  procedure  on  data  from  the  North  Atlantic. 
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under  grant  N00014-86-K-0034. 
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Figure  Captions 


Figure  1.  N-S  density  gradient  vs.  E*W  density  gradient  at  37^  N  and  64^  E  fur 
different  depths.  The  measurements  were  taken  on  August  13-15,  1985  as  part  of  the  REX 
experiment.  The  stars  correspond  to  measurements  over  a  depth  from  105m  to  203m.  Open 
triangles  correspond  to  depths  from  205m  to  603m.  Finally,  dotted  circles  correspond  to 
depths  from  605m  to  743m.  A  stright  line  passing  through  the  origin  with  a  slope  of  -.7 
provides  a  good  fit. 

Figures  2.  Same  data  as  in  Figure  1  segregated  by  depth,  (a)  Upper  layer  (105m  to  203m). 
Note  the  large  scatter,  (b)  Intermediate  layer  (205m  to  603m).  (c)  Bottom  layer  (605m  to 
743m) 


105-203  m 
2  m  resolution 


